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Abstract —We study the secrecy capacity of a helper-assisted 
Gaussian wiretap channel with a source, a legitimate receiver, 
an eavesdropper and an external helper, where each terminal 
is equipped with multiple antennas. Determining the secrecy 
capacity in this scenario generally requires solving a nonconvex 
secrecy rate maximization (SRM) problem. To deal with this 
issue, we first reformulate the original SRM problem into a 
sequence of convex subproblems. For the special case of single¬ 
antenna legitimate receiver, we obtain the secrecy capacity via a 
combination of convex optimization and one-dimensional search, 
while for the general case of multi-antenna legitimate receiver, we 
propose an iterative solution. To gain more insight into how the 
secrecy capacity of a helper-assisted Gaussian wiretap channel 
behaves, we examine the achievable secure degrees of freedom 
(s.d.o.f.) and obtain the maximal achievable s.d.o.f. in closed- 
form. We also derive a closed-form solution to the original SRM 
problem which achieves the maximal s.d.o.f.. Numerical results 
are presented to illustrate the efficacy of the proposed schemes. 


I. Introduction 

The area of physical (PHY) layer security has been pio¬ 
neered by Wyner [I], who introduced the wiretap channel and 
quantified security with the maximal achievable secrecy rate 
(also known as the secrecy capacity) at which the legitimate 
receiver can correctly decode the source message, while the 
eavesdropper can retrieve almost no information. Results in 
El further show that for the classical source-destination- 
eavesdropper Gaussian wiretap channel, the secrecy capacity 
is zero when the quality of the legitimate channel is worse 
than that of the eavesdropper’s channel. One way to achieve 
non-zero secrecy rate in the latter case is to introduce external 
helpers which act as cooperative jammers El- By transmitting 
jamming signals the external helpers degrade the eavesdrop¬ 
per’s channel without hurting the legitimate channel, thus 
allowing secret communication even when the eavesdropper’s 
channel has a much better quality. Works along these lines 
include a-Ei which consider one external helper and ||8l- 
na which consider the case of multiple external helpers. More 
complex relaying scenarios are considered in ifTSl - lfTSl where 
the jamming signal is sent in the relaying phase, or in both 
the broadcasting phase and the relaying phase. 
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Although cooperative jamming approaches improve the se¬ 
crecy rate, their advantage comes from optimally designed 
input covariance matrices, which are difficult to obtain due 
to the nonlinear fractional nature of the problem. To ad¬ 
dress this issue, for the single-antenna eavesdropper case, 
nni-iini propose a suboptimal but cost efficient null-space 
jamming scheme that spreads the jamming signal within the 
null-space of the legitimate receiver’s channel, while ta¬ 
il, M design algorithms to get the optimal solution using 
a combination of convex optimization and one-dimensional 
search. For the multi-antenna eavesdropper case, ina, M 
design the jamming signals so that they align into a pre¬ 
specified jamming subspace at the legitimate receiver, while 
spanning the whole received signal space at the eavesdropper. 
This approach allows the legitimate receiver to completely 
remove the interference by projecting the received signal to 
the secrecy subspace, while confounding the eavesdropper. 
Still for the multi-antenna eavesdropper case, the work of 
PI provides a closed-form expression for the structure of the 
jamming signal covariance matrix that guarantees secrecy rate 
larger than the secrecy capacity of the wiretap channel with 
no jamming signals. The results of HI are obtained under the 
power covariance constraint. 

In this paper, we consider a multi-input multi-output 
(MIMO) Gaussian wiretap channel with one external multi¬ 
antenna helper as in I?]. Different from H, we investigate the 
secrecy rate maximization (SRM) problem under an average 
power constraint. To the best of the authors’ knowledge, 
determining the exact secrecy capacity of a helper-assisted 
MIMO Gaussian wiretap channel has not been previously 
addressed. We first address the problem for the special case 
of single-antenna legitimate receiver. By decomposing the 
original nonconvex SRM problem into a sequence of convex 
subproblems, we obtain the optimal solution to the original 
SRM problem via a combination of convex optimization and 
one-dimensional search. For the general case of multi-antenna 
legitimate receiver, we propose an iterative algorithm to solve 
the original SRM problem via employing the Gauss-Seidel ap¬ 
proach, which successively optimizes each variable while the 
other variables are kept fixed. Specifically, each subproblem is 
convex and admits an optimal solution. Though the proposed 
iterative algorithm provides no guarantee of finding the global 
optimal solution, it constitutes an efficient way for attaining a 
meaningful achievable secrecy rate. 

In order to gain more insight into how the secrecy capacity 
of a helper-assisted MIMO Gaussian wiretap channel behaves, 
we examine the rate at which the secrecy capacity scales with 
log(P), i.e., the maximal achievable secure degrees of freedom 
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(s.d.o.f.) lfT9ll . To this end, we first introduce an alternative 
optimization problem, i.e., maximizing the dimension of the 
subspace spanned by the message signal received at the 
legitimate receiver, under the constraints that the message 
and jamming signals lie in different subspaces at the legiti¬ 
mate receiver, but are aligned into the same subspace at the 
eavesdropper. We then give a critical lemma, proving that the 
maximal achievable objective value of the newly introduced 
optimization problem equals the maximal achievable s.d.o.f.. 
Consequently, the original s.d.o.f. maximization reduces to 
the newly introduced optimization problem. Subsequently, we 
solve analytically the newly introduced optimization problem, 
thus obtaining the maximal achievable s.d.o.f. of the helper- 
assisted MIMO Gaussian wiretap channel in closed-form. 
Further, we derive an analytical solution to the original SRM 
problem, which achieves the maximal s.d.o.f.. Our analytical 
results prove that for the special case of single-antenna legit¬ 
imate receiver, a s.d.o.f. of 1 can be achieved if and only if 
iVe < Na + Nj — 1; for the case of multi-antenna legitimate 
receiver, the maximal achievable s.d.o.f. is zero if and only if 

Ne > Na + Nj. 

We should note that the s.d.o.f. for the helper-assisted 
Gaussian wiretap channel has also been investigated in m, 
II 20 I - II 22 II . Different from our work, the work of Cl studies a 
scenario in which a large number of helpers is available, and 
exploits multiuser diversity via opportunistic helper selection 
to enhance security. The works of GO), im consider a 
scenario in which each terminal is equipped with one antenna, 
while the work of ll^ considers the special scenario in 
which the source, the legitimate receiver and the eavesdropper 
are equipped with the same number of antennas. Further, 
the works of Il20ll - ll22l examine the s.d.o.f. based on real 
interference alignment, while our work is based on spatial 
interference alignment. 

The rest of this paper is organized as follows. In Section II, 
we describe the system model for the MIMO Gaussian wiretap 
channel with one external multi-antenna helper, and formulate 
the secrecy rate maximization problem. In Section III, we 
consider the special case of single-antenna legitimate receiver. 
We investigate the secrecy rate maximization problem, and 
examine the conditions under which a secure degrees of 
freedom equal to 1 can be achieved. In Section IV, we consider 
the general case of multi-antenna legitimate receiver, inves¬ 
tigate the secrecy rate maximization problem, and examine 
the maximal achievable secure degrees of freedom. Numerical 
results are provided in Section V and conclusions are drawn 
in Section VI. 

Notation: A^, tr{A} and rank{A} stand for the hermitian 
transpose, trace and rank of the matrix A, respectively; A(;, j) 
indicates the j-th column of A while and A(;,f : j) denotes 
the columns from i to j of A. span(A) and span(A)'*' are 
the subspace spanned by the columns of A and its orthogonal 
complement, respectively; null(A) denotes the null space of 
A; span(A)/span(B) = {x|x S span(A), x ^ span(B)}. 
A ^ B denotes that A —B is a hermitian positive semidefinite 
matrix. indicates a N x M complex matrix set. i G Z 

denotes that i is a positive integer. I represents an identity 
matrix with appropriate size. Besides, a~^ = max(a, 0); |a| 


is the magnitude of a; a; ^ CAf{0, S) means x is a random 
variable following a complex circular Gaussian distribution 
with mean zero and covariance E. 

II. System Model and Problem Statement 

We consider the MIMO Gaussian wiretap channel with a 
cooperative jammer (see Fig.l) where the source, the legit¬ 
imate receiver, the eavesdropper and the external helper are 
equipped with Na, Nb, Ne and Nj antennas, respectively. The 
source wishes to send its message, x, to the legitimate receiver, 
without being eavesdropped by the eavesdropper. Towards that 
objective, the source is aided by a cooperative terminal, which 
simultaneously transmits jamming signal, z, to confuse the 
eavesdropper. The signals received at the legitimate receiver 
and the eavesdropper can be respectively expressed as 

yd = HiVx + G2WZ-f rid (la) 

Ye = GiVx + H2WZ + rie, (lb) 

where V and W are the precoding matrices at the source and 
the helper, respectively; rid ^ CA/^(0,I) and rie CA/^(0,I) 
represent noise at the legitimate receiver and the eavesdropper, 
respectively; G 2 G and H 2 G represent 

the helper to legitimate receiver and the helper to eaves¬ 
dropper channel matrices, respectively; Hi G 
Gi G denote the channel matrix from the source 

to the legitimate receiver and the source to the eavesdropper, 
respectively. All channels are assumed to be flat fading. 
We assume that global channel state information (CSI) is 
available, including the CSI for the eavesdropper. This is 
possible in situations in which the eavesdropper is normally an 
active member of the network, communicating nonconfidential 
information with the other parties in other time slots a. A 
minimum-Mean-Square-Error (MMSE) receiver is considered 
at the legitimate receiver and the eavesdropper. The rate at the 
legitimate receiver and the eavesdropper can be respectively 
expressed as 

i?d = log|I -k (I -k G 2 QiGf )-iHiQ,Hf I (2a) 

Re = log|I + (I + H2QiHf )”iGiQ„Gf I, (2b) 

where Qa = VV'^ and = WW'^ are the transmit 
covariance matrices at the source and the helper, respectively. 
In the paper, we focus on the SRM problem ll2^ . i.e., Q 

Cs = max i?d — Re 
{QahO,QjhO} 

s.t. tr{Qa -k Qj} < P, (3) 

where P is a given total transmit power budget and Cg denotes 
the maximal achievable secrecy rate, also known as the secrecy 
capacity. 

Generally, the optimization problem of (3) is nonconvex. 
It is challenging and still an open problem to determine the 
exact secrecy capacity. In this paper, we propose to solve the 
problem of (3) by reformulating it into a sequence of convex 
problems. Also, we study the rate at which the secrecy capacity 

*For a given point {Qa,Qj}, the achieved secrecy rate is max(i?^ — 
i?e,0). For ease of exposition, the trivial case with zero achievable secrecy 
rate is omitted. 
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Fig. 1: MIMO wiretap channel with an external helper 


scales with log(P), i.e., the maximal achievable s.d.o.f. lfT9l . 
which equals 


s.d.o.f = 


lim 

P—^OO 


log P' 


(4) 


We compute s.d.o.f analytically and determine its connection 
to system parameters, i.e., the number of antennas at each 
terminal. 

In the following sections, we will begin with the simple case 
where the legitimate receiver is equipped with one antenna. 
Then, a more complicated scenario, in which each terminal is 
equipped with multiple antennas will be investigated. 


III. Helper-Assisted MISOME Wiretap Channel 
In this section, we consider the helper-assisted multi-input 
single-output multi-antenna-eavesdropper (MISOME) wiretap 
channel where the legitimate receiver is equipped with a single 
antenna (Nf, = 1). In such case, the legitimate receiver can 
receive at most one data stream. Thus, the source transmits 
one data stream x. Let v denote the precoding vector at the 
source. The signals received at the legitimate receiver and the 
eavesdropper can be respectively expressed as 

yd = hivx + g2Wz + Ud (5) 

ye = Givt -f H 2 WZ -f rie, (6) 

where hi £ denotes the channel vector from the source 

to the legitimate receiver, and g 2 £ represents the 

channel vector from the helper to the legitimate receiver. The 
rate at the legitimate receiver and the eavesdropper can be 
simplified as, 

Rd = log(l + (IH- g 2 Qjgf )“^hivv^hf) (7a) 

Re = log(l + v^Gf (I + H 2 Q,Hf )"'Giv). (7b) 

Correspondingly, the secrecy capacity equals 
C's = max Rd — Re 

{v.Qi^O} 

s.t. tr{vv^ + Qj} < P- (8) 

Due to the presence of the multi-antenna eavesdropper, 
the SRM problem becomes more complex as compared with 


the problem considered in a-a, cni. To cope with this 
issue, we resort to theorems on partial ordering of hermi- 
tian matrix ll24l . Based on these theorems, we transform 
the original matrix inverse constraint into a convex linear 
matrix inequality (EMI) constraint, which, together with the 
semidefinite relaxation (SDR) technique, enables us to recast 
the original nonconvex optimization problem into a sequence 
of semidefinite programmes (SDPs). Eurther, we prove that 
the optimal solutions to the relaxed optimization problem 
are also optimal solutions to the original SRM problem. 
Consequently, we obtain the optimal solution to the original 
SRM problem with a combination of convex optimization 
and one-dimensional search. On the other hand, to gain more 
insight into how the secrecy capacity of the helper-assisted 
MISOME Gaussian wiretap channel behaves, we examine the 
conditions under which a s.d.o.f. equal to 1 can be achieved. 
To this end, we first introduce an alternative optimization 
problem which keeps the message signal and the jamming 
signal into different subspaces at the legitimate receiver, but 
aligns them into the same subspace at the eavesdropper. We 
then give two key lemmas. Lemma 1 proves that, a s.d.o.f. 
equal to 1 can be achieved if and only if the newly introduced 
optimization problem returns a nonempty set. Lemma 2 gives 
the conditions under which the newly introduced optimization 
problem returns a nonempty set. Combining the two lemmas, 
we finally show that a s.d.o.f. equal to 1 can be achieved if 
and only if Ng < Na + Nj — 1. 

A. Secrecy rate maximization 

To solve the SRM problem in (8), the Two-Layer idea of 111 
is adopted. The key insight is to recast the original optimiza¬ 
tion problem in (8) as a two-level optimization problem. The 
inner-level part is dealt with the SDR technique, and the outer- 
level part is handled by one-dimensional search. Specifically, 
the outer-level part is 

max log(l -f 5(r)) - log(l -f r), (9) 

where p(t) is obtained by solving the following inner-level 
part optimization problem for a fixed t; 

5 (t) = max hivv^hf/(I-f ) (10a) 

s.t. v^Gf(H-H 2 QjHf)“^Giv< r (10b) 
-t- Qj} < P. (10c) 

By performing one-dimensional search on t, the optimal 
T* maximizing the objective function in (9) can be found. 
Correspondingly, the optimal solution {v*, Q*} to the original 
optimization problem of (8) can be obtained. 

In (9), Til, and Tub denote the lower and upper bound on 
7e, respectively. Eirstly, it is obvious that 7e is no less than 
0. Thus, we have Tib = 0. Secondly, according to the security 
requirement, 7e should be no more than 7^. Eurther, 7^ is 
upper bounded by the maximal received signal-to-noise ratio 
(SNR) value of P|hip at the legitimate receiver. Therefore, 

Tub = ^’Ihip- 

So far, Tib and Tub have been determined. In the following, 
we focus on solving the optimization problem of (10), which 
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is still nonconvex. To solve it, we resort to the SDR technique 
of ||25]| . On denoting Qa = vv^ and dropping the rank-one 
constraint, the optimization problem of (10) can be rewritten 
as 

/(r) = max biQ^hf/(I + gaQjgf) (11a) 

{QatO,QjtO} 

s.t. GiQaGf ^ t(I + H 2 Q,Hf) (1 lb) 

tr{Qa + Qj} < P, (11c) 

where the replacement of the constraint (10b) with (11b) can 
be proven using basic theorems on partial ordering ll24l as 
follows: 

(10b) 4^A„,ax((I + HsQ.Hf )”'Givv^Gf) < r 
<^(1 -f HaQjHf )-iGiVV^Gf ^ rl 
<^Givv^Gf ^ t(I -t- HaQjHf) (11b). 

In the above, Amax(A) denotes the maximum eigenvalue of 

A. 

Letting C = (1 + g 2 Qjgf > 0, Qa = CQa, 
and using the Charnes-Cooper transformation 12^ . we can 
recast the optimization problem of (11) as 

/(r) = max hiQahf 

{QatO.QjtO.J>0} 

s.t. ^ + g 2 Qigf = l 

GiQaGf ^T(CI + H2Q,Hf) 

tr{Qa + Qj} < ^P, (12) 

which is a SDP and can be efficiently solved using available 
software packages, e.g., CVX ll26l . 

Let us consider the power minimization problem associated 
with (11), which can be formulated as follows: 

min tr{Qa + Q,} 

{Qa,Qi} 

s.t. hiQahf /(1 + g 2 Qjgf) > /(r) 
GiQaGf ^T(I + H2Q,Hf) 
Qa^O,Q,hO, (13) 

where /(r) is obtained by solving the optimization problem 
of (12). We have the following two propositions. 

Proposition 1: Denote the optimal solution to (13) as 
{Qa,Qj}. Then, Qa is rank-one provided that a positive 
secrecy rate is achieved. 

Proof: See Appendix lAl ■ 

Proposition 2: Denote the optimal solution to (13) as 
{Qa,Qj}- Then, {Qa,Qj} is also the optimal solution to 
the problem of (11). 

Proof: See Appendix iBl ■ 

Let Q° = Qa and Q° = Q j. Combining Proposition 1 with 
Proposition 2, we get that {Q®, Q^} is the optimal solution to 
the problem of (11), such that rank{Q°} = 1. Therefore, the 
optimization problem of (11) is indeed a tight approximation 
of the optimization problem of (10). Moreover, {Q°,Q°} is 
also the optimal solution to the problem of (10) and p(t) = 

/W- 


B. Conditions to ensure s.d.of. equal to 1 

As stated in the preceding sections, it is difficult to obtain an 
analytical expression for the secrecy capacity for the helper- 
assisted Gaussian wiretap channel. Instead, in this subsection, 
we investigate the conditions under which a s.d.o.f. equal to 1 
can be achieved. To this end, we first introduce an alternative 
optimization problem as follows: 


{v,W} 

(14a) 

span(Giv) C span(H 2 W) 

(14b) 

span(g 2 W) n span(hiv) = {0} 

(14c) 

hiv > 0. 

(14d) 


Specifically, we aim to find the point at which the subspace 
spanned by the message signal and that spanned by the 
jamming signal have no intersection at the legitimate receiver, 
such that Rd scales with log(P). Simultaneously, the sub¬ 
space spanned by the message signal belongs to the subspace 
spanned by the jamming signal at the eavesdropper, such that 
Re converges to a constant as P approaches to infinity. 

In the sequel, we first give two key lemmas. Lemma 1 
proves that s.d.o.f. equal to 1 can be achieved if and only if the 
optimization problem of (14) returns a nonempty set. Lemma 
2 gives the conditions under which the optimization problem 
of (14) returns a nonempty set. Combining the two lemmas, 
we finally obtain the conditions to ensure s.d.o.f equal to 1 in 
the helper-assisted MISOME Gaussian wiretap channel. 

Lemma 1: The secure degrees of freedom equal to 1 can 
be achieved if and only if the optimization problem of (14) 
returns a nonempty set. 

Proof: Clearly, if the optimization problem of (14) returns 
a nonempty set, then s.d.o.f. equal to 1 can be achieved. So 
the sufficiency holds true. 

We now turn to prove the necessity by contradiction. If the 
optimization problem of (14) returns an empty set, then at 
least one of the constraints in (14) does not hold true. We test 
(14b)-(14d) one by one: 

1) If (14b) does not hold true, then there exists a direction 
along which the eavesdropper can extract the message 
signal without interference, so the rate at which Re 
scales with log(P) is 1. Together with (4),(8) and the 
fact that the rate at which Rd scales with log(P) is at 
most 1 for the multi-input single-output (MISO) source- 
receiver channel, we arrive at s.d.o.f. — 0. 

2) If (14c) does not hold true, then the message signal is 
aligned in the subspace spanned by the jamming signal, 
so Rd converges to a constant when P approaches to 
infinity, which indicates that s.d.o.f. = 0. 

3) If (14d) does not hold true, then |hiv| = 0, which 
indicates that Rd = 0, thus s.d.o.f. = 0. 

Summarizing, if the optimization problem of (14) returns an 
empty set, s.d.o.f. = 0. Therefore, if s.d.o.f. = 1, the 
optimization problem of (14) returns a nonempty set. This 
completes the proof. ■ 

Before proceeding to Lemma 2, we hrst introduce the 
generalized singular value decomposition (GSVD) transform, 
which provides the basis for the proof of Lemma 2 to follow. 


5 


Definition 1 (GSVD Transform): Given two matrices H £ 


CtVxM Q g ^NxK^ jgj 


k =rank{[H^,G^]'^} 

(15a) 

p =dim{span(H)^ fl span(G)} 

(15b) 

r =dim{span(H) fl span(G)^} 

(15c) 

s =dim{span(H) fl span(G)}, 

(15d) 

then we have 


k = min{M + K, A} 

(16a) 

p =k — min{M, A} 

(16b) 

r =k — min{ A, A} 

(16c) 

s = min{M, A} -f min{ A, A} — k. 

(16d) 

The proof is given in Appendix O According to 
GSVD of (H^,G^) returns unitary matrices 'J'l 

na, the 

^ ^MxM 


and ’S'2 £ i^kxk^ non-negative diagonal matrices Di £ 
£Mxk g f^Kxk^ ^ matrix X £ with 

rank{X} = k, such that 

MS'! = XDf (17a) 

G^2 = XDf, (17b) 



Ir 

0 

0 ■ 


■ 0 

0 

0 

in which Di = 

0 

Si 

0 

, D2 = 

0 

S2 

0 


0 

0 

0 


0 

0 

Ip. 


where the diagonal entries of Si £ and S2 £ are 

greater than 0, and Df Di -|- Df D2 = I. 


For simplicity, in the following text of this paper, we denote 
the above GSVD Transform as 

(^1, ’S'2, Di, D2, X, k, r, s,p) = gsvd(H^, G^). 

Lemma 2: The optimization problem of (14) returns a 
nonempty set if and only if Ng < Na + Nj — 1. 

Proof: We start with the constraint of (14c). With the 
GSVD Transform of (hf, gf), we get si = dim{span(hi) fl 
span(g2)} = 1. Thus to satisfy (14c), we should have |hiv| = 
0 or |g2W| = 0. However, |hiv| = 0 contradicts (14d), so 
we should have |g2W| = 0. 

Without loss of generality, let W = FWi where F = 
null{g2} £ Wi £ C(^j-i)x(tvi-i), Substitut¬ 

ing W = FWi into (14b), we arrive at 

span(Giv) C span(H2FWi), (18) 

in which Gi £ C^-xjv„ a H2r £ cN,xiN^-i)_ 

Invoking the GSVD Transform of (Hf, Gf), we obtain 

(’S'l, ’S'2,Di,D2,X,fc2,'r2,S2,P2) = gsvd(Hf, Gf), 

such that 

H2 4'i = XDf (19a) 

Gi’S '2 = XDf, (19b) 

where ^2 = minjiVa + Nj — 1, Ng}, p 2 = ^2 ~ minliVj — 
r 2 = k 2 - uim{Na,Ng} and S2 = uim{Na,Ng} + 
minjiVj — 1, Ng} — ^2- 


To satisfy (14d), we should have v 7^ 0, which, together 
with (19a) and (19b), indicates that (18) holds true if and 
only if p 2 <Na- 

1) For the case of Ng < Nj — 1, p 2 = Ng — Ng = 0. So 
P2 < Na- 

2) For the case of Nj — I < Ng < Na + Nj — 1, p 2 = 
Ng - (Nj - 1) < Na. 

3) For the case of Ng > Na + Nj — I, P 2 = Na + Nj — 
1 - (Nj - 1) = Na. 

Summarizing, p 2 < Na if and only if Ng < Na + Nj — 1. 

Therefore, if Ng < Na -f Nj — 1, (18) holds true, thus 
the optimization problem of (14) returns a nonempty set. 
Otherwise, if Ng > Na + Nj — 1, (18) does not hold true, 
so the optimization problem of (14) returns an empty set. In 
a word, the optimization problem of (14) returns a nonempty 
set if and only if Ng < Na + Nj — 1. This completes the proof. 

■ 

Throrem 1: The secure degrees of freedom equal to 1 can 
be achieved if and only if Ng < Na -\- Nj — 1. 

Proof: Combining Lemma 1 with Lemma 2, it is clear 
that the secure degrees of freedom equal to 1 can be achieved 
if and only if Ng < Na + Nj — 1. This completes the proof. 


Remark: It is worthwhile to note that Lemma 2 also provides 
us with a way to determine the precoding matrices at the 
source and the helper to achieve s.d.o.f. of 1. In the remaining 
text of this subsection, we first give the precoding matrices 
to achieve s.d.o.f. of 1 in closed-form. We then substitute the 
derived precoding matrices into (8) and solve for the optimal 
power allocation between the message signal and the jamming 
signal. Consequently, we obtain an analytical lower bound on 
the secrecy capacity. 

Revisiting the proof of Lemma 2, for the case of Ng < 
Na + Nj - 1, ^2 = Ng, r 2 = k 2 - min{7Vo, Ng} and S 2 = 
min{7Vo, 7Ve}-|-min{7Vj — 1, A’"e} —^2- Actually, S 2 > 0 which 
can be justified as follows: (i) For the case of Ng < Nj — 1, 
S 2 = min{Wa, Ng} > 0; (ii) For the case of — 1 < A^e < 
Na + Nj -1, S 2 = min{A„ Ng} - (Nj - 1) - Ng > 0 
min{ Aa, Ng} > Ng — (Nj — 1) where the latter inequality can 
be easily verified. 

Let I = {j\r 2 + 1 < j < r 2 + S 2 ,_j £ Z}._Since S 2 > 
0, I is a nonempty set. Let WJ = i)| and 

Vo = ^2(:,* + Na - /C2)/|^2(:,* + Na - k 2 )\ where i £ I. 
According to (19a) and (19b), we arrive at span(GiVo) = 
span(H2FW°) = span(X(:, j)), so {vo,FW°} is a feasible 
point to (14). Substituting v = y/xVo and W = \/P — a;FW° 
into (8), 0 < X < P, we arrive at 


a 


-sub 


l-f|hiVopa; 

max log ;-r 

{0<x<P} ^ 1 -f 


( 20 ) 


in which 


7 


sub 

e 


= Gf 

= a;|GiVo| 
= a;|GiVo| 


(I + (P - x)H2W°W“^Hf )■ 
2 a;(P-a;)|W?"HfGiVof 
1-f (P-a;)|H2W°|2 

2 a:(P-a;)|H2W?|2|GiVof 


Givo 

(21a) 


1 + (P-x)|H2WJ|2 


(21b) 
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where (21a) follows from the matrix inverse lemma m, 
and (21b) follows from the fact that span(GiVo) = 
span(H 2 W°) = span(X(:, i)). 

For ease of exposition, let a = jhiVop, b = |GiVop and 
c = IH 2 WJP. Also, noting that the logarithm function is a 
monotonically increasing function, therefore the optimization 
problem of (20) becomes 

= max rjix), (22) 

0<x<P 


in which 


1 + ax 

1 + bx — [bc{P — x)xl (1 + c(P 
(1 + aa;)[l + c{P - a:)] 

1 + cP + (6 — c)x 


a^))] 


(23) 


Resorting to carefully mathematical deductions, we solve (22) 
and arrive at that when the total transmit power P is big 
enough. 


Cr” ~log(aP) - 21og(l + y/V^), (24) 


Applying Lemma 3, we arrive at, respectively. 


In |(I + G 2 QjG^) = max (^t,(So) (27a) 


In 


(I + H2Q,Hf+ GiQaGf) 


maxi^e(Si), (27b) 


where (^t,(So) = -tr{So(I +G 2 QjG^)} + In |So| + iVf,, and 

(^e(Si) = -tr{Si(I+H 2 QjHf+GiQ,Gf)}+ln|Si|+7Ve. 

Substituting (27a) and (27b) into (2a) and (2b), respectively, 
we arrive at 


Rd = max(/Jb(So) + In |l + HiQaHf + G 2 QjGf^| (28a) 

Sq^O 

Rf. = — max (^e(Si) — In |l + H 2 QjH^| . (28b) 

Si ^0 

Further, substituting (28a)(28b) into (3), we arrive at 


Cs 

s.t. 


max 

{Qa^O,Q3^0,So^O,Si^O} 


0 (So,Si,Q,,Q,) 


tr{Qa + Qj} < P, 


(29) 


where the details are given in Appendix D. 

Substituting a = jhiVop, b = |GiVop and c = |H 2 W°p 
into (24) yields 

Cr” ~log(|hiv„|2p) - 21og(l + ^|Giv„| 2 /|H 2 W?P), 

(25) 

which in turn explicitly corroborates that s.d.o.f. equal to 1 
has been achieved. 

IV. Helper-Assisted MIMOME Wiretap Channel 

In this section, we consider the helper-assisted MIMOME 
wiretap channel where each terminal is equipped with multiple 
antennas. In such MIMO case, the SRM problem becomes 
more complex as compared with the problem considered in 
Section III, and actually, it is still an open problem. To deal 
with this issue, we hrst reformulate the SRM problem in (3) to 
a form that can be processed with the Gauss-Seidel approach, 
which successively optimizes each variable given that the 
other variables are fixed, thus giving an iterative algorithm 
to solve (3). We then examine the maximal achievable s.d.o.f. 
and reveal its connection to the number of antennas at each 
terminal. We obtain both the maximal achievable s.d.o.f. and 
the solution that achieves the maximal s.d.o.f. in closed-form. 

A. Secrecy rate maximization 

In order to reformulate the SRM problem in (3) to a form 
that can be processed with the Gauss-Seidel approach, we need 
the following lemma. 

Lemma 3: Given a positive definite matrix E G it 

holds that 

ln|E“^|= max (pfS), (26) 

se'^x'^.s^o 

where (/?(S) = —tr(SE)-t-ln |S|-|-W Moreover, for the optimal 
solution to the right-hand side of (26), it holds that S* = E~^. 

Proof: Please refer to ||28]. ■ 


where ^(So, Si, Qa, Q^) = (/^(So)-f V5e(Si)-t-w(Qa, Qj) 
in which w(Qa,Qj) = In |l-f HiQaHf^-f G 2 QjG|^| -f 
Injl-t-HaQjHf |. 

Although the optimization problem of (29) is still not con¬ 
vex, we observe that if we fix either {Qa, Qj} or Sfii = 1, 2), 
the remaining problem is convex and can thus be solved 
efficiently. Hence, we turn to a two-stage iterative method 
(Gauss-Seidel approach), and approximately solve the opti¬ 
mization problem of (29) via iterations between the following 
two subproblems. 

1 ) Eix {Qa,Qj}, and maximize 0(So, Si, Qa, Q^) over 
{So,Si}. 

2) Eix {So, Si}, and maximize 0(So, Si, Qa, Qj) over 

{Qa,Qj}. 

Specifically, when (Qa, Qj} is fixed, let 

{So,S{} = arg max 6»(So, Si, Qa, Qj). 

{So,Si} 

Applying Lemma 3, we arrive at 

SS = (I + G2QiGf)-i, 

S} = (I-f H2QjHf + GiQaGf 

Besides, when {So,Si} is fixed, the maximization of 
6 »(So,Si,Qa,Qj) over {Qa,Qj} is a convex optimization 
problem and can be efficiently solved by available software 
packages, e.g., CVX 1261. 

One can easily verify that the above iterative process leads 
to a monotonically non-descending objective function value 
of 0(So, Si, Qa, Qj). Moreover, for a given limited transmit 
power, the achievable secrecy rate is upper bounded. Thus, the 
above iterative algorithm is convergent. 

Remark: Although the above iterative algorithm provides no 
guarantee of finding the global optimal solution to the problem 
of (29), our numerical results in the following section show 
that it attains a fairly good secrecy rate performance. 
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B. Maximal achievable secure degrees of freedom 

In this subsection, we examine the maximal achievable 
s.d.o.f. and determine its connection to the number of antennas 
at each terminal. Similar to Section III, we first introduce an 
alternative optimization problem as follows; 

d = max rankjHiV} (30a) 

{v,w} 

s.t. span(GiV) C span(H2W) (30b) 

span(G2W) n span(HiV) = {0}. (30c) 

Specifically, we find the feasible points at which the subspace 
spanned by the message signal and that spanned by the 
jamming signal have no intersection at the legitimate receiver. 
Simultaneously, the subspace spanned by the message signal 
belongs to the subspace spanned by the jamming signal at the 
eavesdropper. Among these feasible points, we determine the 
one at which the rank of HiV is maximized. 

Lemma 4: The maximal achievable secure degrees of free¬ 
dom, defined in (4), equal to d. That is, s.d.o.f = d. 

Proof: See Appendix |E] ■ 

From Lemma 4, we observe that to obtain the s.d.o.f, we 
need only to focus on solving (30). To this end, in the sequel, 
we first give a heuristic method which gives a closed-form 
feasible point {V,W} to (30), followed by the derivation of 
d* = rank{HiV} in closed-form. Subsequently, in Lemma 5, 
we prove that d = d*. Combining Lemma 4 and Lemma 5, we 
finally obtain s.d.o.f = d*. Further, we prove that {V, W} 
achieves the maximal s.d.o.f., i.e., {V,W} constituting the 
s.d.o.f.-optimal solution to the original SRM problem. 

The aforementioned heuristic method to obtain {V, W} is 
shown in Table 1. Notice that in Table I, null{Gi} returns an 
empty matrix when Na < Ne. In the following text, we prove 
that {V,W} is a feasible solution for (30), and derive the 
closed-form expression for d*. As in Table I, four cases are 
discussed. 

In Case I and Case II, it is clear that {V, W} is feasible 
to (30) and d* = rank{HiV} = minjA’a, Nt}. 

In Case III, for the subcase of do + di > Nb, V = 
[Vo,Vi] and W = Wi. According to (31), we get 
span(G2Wi) = {0} and span(H2Wi) = span(GiVi). 
In addition, GiVq = 0. So, span(H2W) = span(GiV) 
and span(HiV) Pi span(G2W) = {0}, which indicate that 
{V,W} is feasible to (30). Furthermore, Vq is orthogonal 
with Vi by definition, thus 

d*=rank{[Vo,Vi]} 

=rank{Vo} -f rank{Vi} = Nb. 

For the subcase of dp-l-di < Nb, V = [Vq, Vi, V2] and W = 
[Wi, W2]. As in the subcase of do + di > Nb, GiVq = 0 , 
span(G2Wi) = {0} and span(H2Wi) = span(GiVi). In 
addition, according to (32), span(H2W2) = span(GiV2) 
and d2 = min{s4, Thus span(H2W) = 

span(GiV) and span(HiV) P span(G2W) = {0}. There¬ 
fore {V,W} is feasible to (30). Furthermore, [Vo,Vi] is 


TABLE I; A heuristic method to obtain {V,W} which is 
_feasible to (30) _ 

Case I: Na > Ne + Nfj. Let V = null{Gi} and W = 0 . 

Case 11 : Nj > N + Ne. Let W = null{G2} G c^aX(Nj-Nk) 

V be the right singular matrix of Hi. 

Case III: Na < Ne + A^b and A^b < < A^e + A^b- For a start, let 

Vq = null{Gi} and do = {Na — A^e)"*"- Secondly, denote H2 = 
H2r g_C^ex(JV 7 -^ 6 ) where T = null{G2} G 
Denote Gi = GiVq where Vg = nulljVQ } G Invoking 

the GSVD Transform of (H|^, G^) yields 

(^i,^ 2 ,Di,D 2 ,X,fe 3 ,r 3 ,S 3 ) = gsvd(Hf ,Gf). (31) 

Subsequently, let di = 53, C3 = ra + Ne — ^3, and check 

1 ) If do + c/i > Nf,, let Wi = (:, 7*3 + 1 : 7*3 + Nf, — do) and 

Vi = Vg^2(:, C3 + I : C3 + Nb-do). Lastly, let V = [Vq, Vi] 
and W = Wi. 

2 ) Otherwise, let Wi = r^'i(:,r3 + 1 : 7^3 + S3) and 

Vi = Vg^2(:,C3 + 1 : C3 + S3). Thirdly, denote = 

null{[Vo, G c^“X(iV„-tio-di) and Gi = GiVg^. 

Invoking the GSVD Transform of (H|^, Gj^) yields 

(^i,4'2,Di,D2,X,fc4,r4,M) = gsvd{Hf,Gf). (32) 

Then let W2 = 4'i(:,r4 -|- 1 : 7-4 -|- ^2) and V2 = Vqj^^2(: 

, C4 + 1 : C4 + d^) in which d2 = min{s4, [— — (^o+'Ii) j| and 
C4 = r4 -I- (Na — do — di) — k4. Lastly, let V = [Vq, Vi, V2] 
and W = [Wi, W2]. 

Case IV: Na < N^-i-Nb and Nj < Nb- For a start, let Vq = null{Gi} 

and do = (Na — Ne)'^. Secondly, denote Vg = null{Vg^} G 

and Gi = GiVg. Invoking the GSVD Transform of (Hj^jCj^) yields 

(^i,^2,Di,D2,X,fc4,r4,S4) =gsvd(H|^,Gf). (33) 

Then let W2 = ^’i(:, -\- 1 ■. -\- ^2) and V2 = Vq^'2(:, C4 -|- 1 : 

C4 + df) in which d^ = min{s4, g ^ *^0 J j. and C4 = r4 -f Na — k^. 
Lastly, let V = [Vq, V2] and W = W2. 


orthogonal with V2 by definition, thus 

d’^=rank{[Vo,Vi,V2]} 

=rank{[Vo, Vi]} -|-rank{V2} 

=rank{Vo} -f rank{Vi} -f rank{V2} 

= min{do -\- di-\- d 2 , Na}. 

In Case IV, V = [Vo,V2] and W = W2. Accord¬ 
ing to (33), span(H2W2) = span(GiV2), which, together 
with GiVq = 0, gives span(H2W) = span(GiV). In 
addition, span(HiV) fl span(G2W) = {0} due to d 2 = 
min{s4, Thus, {V,W} is feasible to (30). Fur¬ 

thermore, Vq is orthogonal with V2 by definition, therefore 

d* =rank{[Vo,V2]} 

=rank{Vo} -f rank{V2} 

= min{do -f d 2 ,Na}. 

Summarizing the above four cases, we can rewrite d* into 
a more compact form as follows: 

d* = min{do -f d* -f d^, Na, Nt], 


(34) 








TABLE II: Summary of the closed-form results on d*(s.d.o. f.) 


Inequalities on the number of antennas at terminals 

d* (s.d.o.f.) 

Na>Ng+ Nb 

m'm{Na,Nb} 

Nj >Ng + Nb 

2Nb + Ng- Nj <Na<Ng + Nb 

Nb < Nj <Ng + Nb 

Nb + Ng- Nj <Na< 2Nb -yNg- Nj 

Nb < Nj <Ngy-Nb 

Na + Nj - (Nb + Ng) + min{s, ^2iV,+iVe-iV„-iV, 
s = min{Nb Ng — Nj,Ng} -|- mm{Nj,Ng} — Ng 

Ng<Na<Ng + Nb, Nj < Nb 

Na Ng-y minis, J “ min{Arj,Are} 

Na<Nb + Ng- Nj, Nb < Nj <Ng + Nb 

minis, L^J} 

s = min{Afa, A^e } -I- min{At,- ,Ng} — min{A^a + Nj, Ng } 

Na < Ng, Nj < Nb 


in which 

d*o = {N, - Ne)+ (35a) 

dl = imm{Na, N,} + {N, - Nb)+ - Ne)+ (35b) 

d* = min{s, ([— -)+}^ (35c) 

where s = min{A^a — (do + dl),Ne} + min{A(j, A"e} — 
min{Ara - (dj + d^) + A^-, A^e}- 

To gain more insight into d*, we give Table II which 
clarihes the connection of d* to the number of antennas at 
each terminal. 

Lemma 5: On d defined in (30), we have d = d* where d* 
is given in (34). 

Proof: See Appendix |F] ■ 

Remark: According to Lemma 5, it is straight-forward that 
the feasible solution {V, W} given in Table I is also the 
optimal solution to (30). 

Throrem 2: Consider a helper-assisted MIMO Gaussian 
wiretap channel, as depicted in Fig.l, where the source, the 
legitimate receiver, the eavesdropper and an external helper 
are equipped with Na, Nb, Ng and Nj antennas, respectively. 
The maximal achievable secure degrees of freedom 

s.d.o.f. = d*, (36) 

where d* is given in (34). 

Proof: Combining Lemma 4 and Lemma 5, it is clear that 
s.d.o.f. = d*. This completes the proof. ■ 

Corollary 1: The feasible point {V, W} for the optimiza¬ 
tion problem of (30), given in Table I, serves as a s.d.o.f.- 
optimal solution to the original SRM problem in (3). It 
achieves the maximal s.d.o.f.. Moreover, Table II clarifies 
the maximal achievable s.d.o.f. of a helper-assisted MIMO 
Gaussian wiretap channel, and reveals its specihc connection 
to the number of antennas at each terminal. 

Proof: With Theorem 2, it is straight-forward to arrive at 
these conclusions. ■ 

Corollary 2: When Nb > 1, the maximal achievable s.d.o.f. 
of a helper-assisted MIMO Gaussian wiretap channel is zero 
if and only if Ng > Na -f Nj. When Af, = 1, the maximal 
achievable s.d.o.f. of a helper-assisted MIMO Gaussian wire¬ 
tap channel is zero if and only if Ng > Na -\- Nj — 1. 

Proof: See Appendix iGl ■ 



Fig. 2; Average secrecy rate versus SNR, Na = Ng = 3, 
Nj =2, Nb = I 

V. Numerical Results 

In this section, we give numerical results to show the secrecy 
rate performance of the proposed schemes and validate our 
theoretical hndings. All channels are assumed to be quasi¬ 
static flat Rayleigh fading and independent of each other, 
with entries distributed as CA(0,1). The noise vector at each 
receiver is assumed to be AWGN, with i.i.d. entries distributed 
as CA(0,1). In each figure, details on the number of antennas 
at each terminal will be depicted. Results are averaged over 
1000 independent channel trials. 

We first test the secrecy rate performance of our proposed 
schemes and compare them with the existing method. 

Fig.2 illustrates the secrecy rate performance of the single¬ 
antenna legitimate receiver case. The lines labeled as “Optimal 
Scheme” and “Alignment Scheme” illustrate the secrecy rate 
performance of Cg and in (8) and (20), respectively. 

The line labeled as “ZF Scheme” [lOl gives the secrecy 
rate achieved by the scheme which completely nulls out the 
jamming signal at the legitimate receiver and matches the 
message signal with the source-receiver (legitimate) channel. 
It shows that both Cg and increase linearly with SNR. 
In contrast, there exists a performance ceiling on the secrecy 
rate achieved by the ZF Scheme. 

Fig. 3 illustrates the secrecy rate performance of the multi¬ 
antenna legitimate receiver case. The bar labeled as “Align- 
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ment Scheme” shows the secrecy rate result of our proposed 
heuristic method. In such case, closed-form precoding matrices 
{V,W} are given in Table I and the total power is equally 
distributed over all message signal streams and jamming 
signal streams. The bar labeled as “Gauss-Seidel Approach” 
shows the secrecy rate performance of our proposed iterative 
algorithm in Section IV. A. As stated in Corollary 1, {V, W} 
given in Table I, actually acts as a s.d.o.f.-optimal solution to 
the original SRM problem in (3). So the initial point is set 
as the closed-form solutions in the Alignment Scheme. For 
comparison, Fig. 3 also plots the secrecy rate performance 
of the method proposed in a, wherein the secrecy rate 
maximization method is derived under a power covariance 
constraint. Thus, to find the maximal achievable secrecy rate 
under an average power constraint, we have to solve [5, 
equation (41)] 1^ 

Rs(P) = max i?s(S). 

S^O,tr{S}<P 

That is, numerical search over the power covariance matrix S 
is performed to compute Rs{P). Since such numerical search 
is based on random choices of S, it is difficult to decide when 
to stop it. To deal with this issue, we first determine the run 
time of our Gauss-Seidel based algorithm of Section IV. A, 
which terminates when the relative secrecy rate improvement 
between two adjacent iterations is less than 10“^. We then 
run the algorithm proposed in Q for the same run time. It 
is worthwhile to note that our proposed Alignment Scheme 
has closed-form solutions, so it is the most computationally 
inexpensive scheme. Fig. 3 shows that, with the same run 
time, our proposed Gauss-Seidel based algorithm achieves 
higher secrecy rate than the algorithm proposed in n. More 
encouragingly, it can be seen that the Alignment Scheme 
achieves nearly the same secrecy rate as the Gauss-Seidel 
Approach based algorithm in the high SNR regime. This is 
consistent with the fact that the Alignment Scheme is s.d.o.f.- 
optimal, thus near-optimal at high SNR. Besides, it can be seen 
that the Alignment Scheme achieves higher secrecy rate than 
the algorithm proposed in Q in the high and medium SNR 
regimes. To gain more insight into the Gauss-Seidel based 
algorithm. Fig. 4 plots the convergence of it. Results show 
that our proposed Gauss-Seidel Approach converges very fast 
and stabilizes after several loops. 

We then test the achievable s.d.o.f. performance for the 
helper-assisted MIMO wiretap channel and validate the theo¬ 
retical results of Section IV. B. 

In Fig. 5, the stem labeled as “Theoretical Results” shows 
the theoretical maximal achievable s.d.o.f. according to Table 
II. The stem labeled as “Numerical Results” shows the s.d.o.f. 
achieved by the proposed Alignment Scheme. In the proposed 
Alignment Scheme, closed-form precoding matrices {V, W} 
are given in Table I and the total power is equally distributed 
over all message signal streams and jamming signal streams. 
The total power P is set as 50dB. For each channel trial, we 
substitute the closed-from solution into (3) and compute the 
secrecy rate C°. We then compute the s.d.o.f. as the rate at 
which the secrecy rate C° scales with log P, i.e., C° /log P. 
It can be seen that the theoretical results almost coincide with 



Fig. 3: Average secrecy rate versus SNR, Na = Ni, = 3, 
Nj =4, Ne = 3 



Fig. 4; Average secrecy rate versus the number of iterations, 
N, =4, Na = Ni, = 3, Ne = 3 



Fig. 5; s.d.o.f. versus Na, Nj = 3, Nb = 3, Ne = 4: 


the numerical results. 

Fig. 6 and Fig. 7 plot the maximal achievable s.d.o.f. for 
the helper-assisted MIMO Gaussian wiretap channel under 
various antenna configurations, according to Table II. Results 
show that for the case of single-antenna legitimate receiver, 
the maximal achievable s.d.o.f. is zero if and only if Ne > 
Na + Nj — 1, while for the case of multi-antenna legitimate 
receiver, the maximal achievable s.d.o.f. is zero if and only if 
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Fig. 6: s.d.o.f. versus Na, Nj = 2, Ne = i 



Fig. 7: s.d.o.f. versus Na, Nj = 4, Ne = 4 


Ne > Na + Nj. Further, it is illustrated that the maximal 
achievable s.d.o.f. benefits from the increasing number of 
antennas at any of the three terminals, i.e., the source, the 
legitimate receiver and the external helper. These results are 
consistent with the theoretical findings of Section IV. B. 

VI. Conclusion 

We have studied the secrecy capacity of the MIMO Gaus¬ 
sian wiretap channel, where a multi-antenna external helper 
is available. For the special case of single-antenna legiti¬ 
mate receiver, we have obtained the secrecy capacity using 
a combination of convex optimization and one-dimensional 
search. For the case of multi-antenna legitimate receiver, we 
have reformulated the original nonconvex SRM problem into 
several convex subproblems. By doing so, we have been able 
to provide an iterative algorithm, which attains a fairly good 
secrecy rate performance. In addition, we have addressed 
the s.d.o.f. maximization analytically. Specifically, we have 
obtained an analytical s.d.o.f.-optimal solution to the original 
SRM problem, based on which, we have obtained the maximal 
achievable s.d.o.f. in closed-form. These results uncovered the 
connection between the maximal achievable s.d.o.f. and the 
system parameters, thus shedding light on how the secrecy 
capacity of a helper-assisted MIMO Gaussian wiretap chan¬ 
nel behaves. Numerical results have validated the theoretical 
findings and confirmed the efficacy of our proposed schemes. 


Appendix A 

Proof of Proposition 1 
The associated Lagrangian of (13) is 

£ =tr{Qa -f Qj} -F /i[/(T)(l -F g2Qigf) - hiQahf ] 

+ tr{Zi[GiQaGf - t(I + H2Q,Hf )]} 

- tr{ZaQa} - tr{ZjQj}, (37) 

where Zi, Z^, Zj and n are dual variables associated with the 
inequalities in (13). The optimization problem of (13) is con¬ 
vex with part of the Karush-Kuhn-Tucker (KKT) conditions 


as follows: 

Za = I-^hfhi-FGfZiGi (38a) 

ZaQa = 0 (38b) 

Zi ^ 0 ,Za ^ 0 ,/i > 0. (38c) 

Substituting (38a) into (38b), we arrive at 

(I-FGfZiGi)Qa = /rhfhiQa. (39) 

Since I -F Gf^ZiGi 0 , so ranklQo} = rank{(I -F 


Gf^ZiGi)Qa} = rank{/ih{^hiQa} < 1. In addition, 
rank{Qa} = 0 implies that = 0, which contradicts the 
positive secrecy rate requirement. Thus, rank{Qa} = 1. This 
completes the proof. 


Appendix B 

Proof of Proposition 2 
For notational simplicity, let 

®(Qa,Qj) = hiQahf/(l-Fg2Qigf). 

Denote the optimal solution to (11) as {Qa, Qj}- Apparently, 
the point {Qa, Qjj is also feasible to (13). Thus we have 

tr{Qa + Qj} < tr{Qa -F Qj} < P, 

which indicates the point {QajQj} is feasible to (11). Ac¬ 
cording to the definition of /(r), we see $(Qa, Qj ) < fir). 
Moreover, from (13), <I>(Qa,Qj) > /(r). Thus, 

^(Qa.Qj) = fir), 

which implies that the point {Qa,Qj} is also the optimal 
solution to (11). This completes the proof. 

Appendix C 

Proof of the equalities in (16) 

Apparently, by the definition of (15a), (16a) holds true. 
According to ll24l . for any given matrix of A, it holds that 

span(A^) = null(A)-L. (40) 

Applying (40) to (15b), we get 

span(H)^ n span(G) = null(H^) n null(G^)^ 

= null(H^)/[null(H^) n null(G")] 

= null(H^)/null([H^,G"]^). (41) 

In addition, nul^p-^, G-^]^) C null(H'^) by definition, so 
we have p = dim{null(H^)} — dim{null([H^, G^]^)} = 
N — min{M, N} — (N — k) = k — min{M, A}. 

Similarly, we can prove (16b)-(16d). This completes the 
proof 
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Appendix D 

Proof of the lower bound on Cs in (24) 


Appendix E 
Proof of Lemma 4 


The optimization problem of (22) can be solved by resorting 
to carefully mathematical deductions. Let y = l+cP+(6—c)x, 
a = min{6, c} and (3 = max(6, c). Then, y G [1 + aP, 1+/3P]. 
Substituting x = (y — (1 + cP))/{h — c) into (23) and making 
some mathematical transformations yield 

Viy) = 1^- {Ay + B/y), (42) 

, ail-\-cP)-\-c — b \2ail-\-cP)-\-c — h]b 

in which n = -^-^, 

c-b [c-bf 

yl = and B = 4i±f£Mi±;a±£j, 

(c -b) (c - b) 

Resorting to (42), the optimization problem of (22) can be 
transformed into a new optimization problem that searches for 
y as follows: 


r]ms.x= max K-{Ay + B/y). (43) 

l+aP<y<l+l 3 P 


1) Lor the case of o(l + cP) > b — c, B > 0. Thus, 
ri{y) is increasing in y when y < yo, and decreas¬ 
ing in y when y > yo. Herein, yo = yjBjA = 
y^5(T+l[P)|a^T+lIP)'^^y+c]7(ac). Therefore, if 1 -f 
OiP < yo ^ 1 + PP, ilmax = K — 2y/AB . Otherwise, 
r]{y) achieves its maximal value at the two endpoints 
(y = 1 -I- aP or y = 1 -f /3P), and p^ax = max{(l -f 
aP)/(l + 6P),l}. 

2) Lor the case of a(l + cP) < b — c, B < 0. Thus, ri{y) 
decreases monotonically with respect to y. It achieves 
its maximal value at the two endpoints, and pmax = 
max{(l -I- oP)/(l -I- bP), 1}. 

Summarizing, if 1 -I- aP < yo < 1 + I3P, then pmax = 
K — 2'/AB and 


cr =iog( 


c(l -I- oP)[c + b — 2\fbc\J 1 -I- (c — 6)/(a -f acP)] 


{c-b) 

c{a — b) + 2(c — a)^/bc^yl + (c — b) / (a + acP) 


{c-b)" 




(44) 


where the optimal solution 

X* = {yo - (1 + cP))/{b - c) 

^ \/6(l + cP)[a(l + cP)-b + c]/{ac) - (1 -f cP) 
b — c 

Otherwise, when a > b, pmax = (1 + ciP)/{l + bP) and 
(^sub _ _|_ aP)/{l + bP), where the optimal solution 

X* = P; when a < b, Pmax = 1 and C™*’ = 0, where the 
optimal solution x* = 0. In the sequel, we refer to these two 
solutions X* = 0 and x* = P as the trivial solutions. 

Consider the nontrivial solution of (45). When the total 
transmit power P is big enough. 


Before proceeding, we first give two critical properties on 
matrix that will be used in the following analyses. That is, for 
any given matrices A and B, if B is invertible, then 

span(A) = span(AB) (47) 

rank{A} = rank{AB}. (48) 

Lirstly, span(A) = span(ABB“^) C span(AB). Secondly, 
span(A) D span(AB). Therefore, the equality (47) holds 
true. With (47), it is clear that the equality (48) holds true. 

Given an arbitrary point of {V, W}, with the definition in 
(4), we can re-express the achieved s.d.o.f. as follows: 

h(V, W) = rank{HiV} - m(V, W) - n(V, W), (49) 

in which m(V,W) = dim{span(GiV)/span(H2W)} and 
n(V, W) = dim{span(G2W) nspan(HiV)}. 

Assume that {V,W} is the optimal solution to (30), then 
we have span(GiV) C span(H2W) and span(G2W) fl 
span(HiV) = {0}. The achieved s.d.o.f. h(V,'W) = 
rank{HiV} = d. In addition, s.d.o.f. > /i(V,W) by 
definition. Thus, s.d.o.f. > d. As such, to complete the proof 
of Lemma 4, we only need to prove s.d.o.f. < d. In the sequel, 
we show that, for any given point of {V, W}, we can always 
find another feasible point for the problem of (30), {V', W'}, 
such that h{V,W) < rank{HiV'} < d, thus giving the proof 
of s.d.o.f. < d. 

Without lose of generality, denote V G (CNaxKa and W G 
With the GSVD Transform of (W^G^,V^Hf), 
we obtain unitary matrices G and 4>2 G 

non-negative diagonal matrices Di G and D2 G 

CKaXks^ and a matrix X G with rank{X} = fcs, 

such that 

G2W’5 'i = XDf (50a) 

HiV 4'2 = XDf, (50b) 

where ko = min{ATo + Kj, Nr}, = k^ — min{ATo, Nr} and 
S5 = min{ Aa, Nr} + min{ Aj, A^.} — k^. 

Let 

#} = 4'i(:,r5-f 1 : rs-f S5) (51a) 

4>o = 1 : rs), 4'i(:,r5 + S5 : A,)] (51b) 

4^2 = 4 '2(:,Ci„-f 1 : Cin-f S5) (51c) 

4 ^° = [ 4 ' 2 (:, 1 : Ci„), 4 ' 2 (:, Qn + S5 + 1 : Ka)], (51d) 

in which Cin = rs -f Ka — k^. Since and 4'2 are invertible 
matrices, 4^^ = and ’S'2 = [^2:’^2] also 

invertible matrices. Applying (47) and (48), we have 

h{Y, W) = h(V%,W%) (52a) 

= rank{HiV4'°}-m(V#f W^'i) (52b) 

<rank{HiV’S'^}-m(V^f W’®';), (52c) 


G™'’«log(aP) - 21og(l-f a/Vc)- (46) ™ '^^ich (52b) can be justified with span(G2W’5''i) n 

span(HiV^ 2 ) = span(HiV^f. Besides, (52c) comes from 
the fact that m(V4'f > m(V#f 


This completes the proof. 
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With the GSVD Transform of 

we obtain unitary matrices G and #2 G 

([^{Ka-s5)x(Ka-s5)^ non-negative diagonal matrices Di G 
(^KjXke y)2 g and a matrix X G 


with rankjX} = k^, such that 

= XDf (53a) 

GiV4'°^2 = XDf, (53b) 

where fcg = minjitia — -\- Kj,Nf.}, tq = — minlitTo — 

S 5 , Ne} and sq = minlitT^, iVe} + minlitTo - 55 , Ng} - fee- 
Let 

= ^i(:,l : re+S 6 ) (54a) 

4'° = «'i(:,r 6 + S 6 (54b) 

= # 2(:,1 : Cin+ S 6 ) (54c) 

4^° = ^2(:, Ci„ + S6 + 1 : Xa - S 5 ), (54d) 

in which Ci„ = re + {Ka — S 5 ) — fee- Since and ’5 '2 are 
invertible matrices, so and ^'2 — 


are also invertible matrices. Applying (47) and (48), we have 

rankjHiV#^} - m(V4'f W#;) 

= rankjHiV’®'^’®'^} - (55a) 

= rank{HiV4'^^^} - rank{^^} (55b) 

< rankjHiV’S'^’CE'^}, (55c) 

where due to span(GiV’S'2’®'y/span(H2W4'j’^2) = 
span(GiV^ 2 ^ 2 )> (55b) holds true. Besides, (55c) holds 
true due to rank{HiV^ 2 ^ 2 } ^ rank{HiV4'f^'2} + 
rankjHiV^^^^} and rank{HiV^^^^} < rankl^f^}. 

Combining (52) with (55), we arrive at 

/i(V,W) < rankjHiV^^^'^}. (56) 

With (53) and (54), we arrive at m(V4'2’^f = 

0, thus 

span(GiV4'°^f C span(H2W4';«'f. (57) 

In addition, with (50) and (51), we get n(V4'f = 0. 

So span(G2W4'2) flspan(HiV’®' 2 ) = {□}, which, together 
with the facts that span(G 2 W’®' 2 ) = span(G 2 W’®' 2 ^j) and 
span(HiV^^) D span(HiV’®'^’^f, gives 

span(G2W^;^;) n span(HiV^°^f = {0}. (58) 

Combining (57) with (58), we know 
is a feasible point for the problem of (30). By definition, 
ranklHiV^o#!} < d, which, together with (56), indicates 
that /i(V,W) < d. 

Because the above derivations hold true for any given point 
of {V,W}, we conclude that s.d.o.f. < d. This completes 
the proof. 

Appendix F 
Proof of Lemma 5 

Clearly, d > d* holds true. So if we can further prove 
d < d*, then the proof of Lemma 5 is completed. In the 


following text, we give the proof of d < d* by contra¬ 
diction. Assume that there exists a feasible point {V,W} 
of (30), where V G d^ = rank{HiV} and 

Ka = d^ > d*. In such case, we have rank{V} = d^ due 
to d'^ = rankjHiV} < rankjV} < Ka = d^. Besides, by 
definition, d^ < min{iVo, Ni,} always holds true. In the sequel, 
we discuss the four cases in Table I and give contradictions 
one by one. 

In Case I and Case II, d* = min{Wo, Nb}. The assumption 
d^ > d* implies d^ > mm{Na, Ni,}, which contradicts the 
fact d^ < minlWo, Nb}. 

In Case III, when d* = minlWo, Nt}, the assumption d^ > 
d* = Nb contradicts the fact d'^ < min{7Vo, Nb}. As such, we 
only need to focus on the case of d* = do + di + d 2 , where 
do = (A^a-./Ve)+, di = S3 and d 2 = min{s 4 , +<^i) j 

1) For the case of 54 < +<^ 1 ) _l_ _l_ 

In addition, d^ > d*. So d** > do + + S4. Therefore 

d** — do > S3 + S4, which contradicts (30b). The 
explanation is as follows. According to (32), S4 = 
min{Aj, TVel+minlWa—do —di, A^e}—min{A(j+A^a— 
do — di, Ng}. With the GSVD Transform of (H|^, G{^), 
s = dim{span(H2) fl span(Gi)} = min{A(,, A^e} + 
min{A"a, Ng] — min{A(j + Na, A^e}- It is easy to verify 
that s < S3 + S4. In addition, to satisfy (30b) we should 
have rank{GiV} < s. Thus rank{GiV} < S3 + S4. 
Moreover, rank{V} — do < min{Ata, Wg} due to the 
fact rank{V} < Na, so rank{GiV} = rank{V} — 
do = d} — do- Therefore, d** — do < S3 + S4, which gives 
the contradiction. 

2) For the case of S 4 > = do + di -f 

|^ Afi,-(do+di) j ^ together with the assumption d* > 

d*, gives 

d# > do + di + + 

If Atf, — (do + di) is an even number, (59) is equivalent to 
2d* > Atft-l-do + di. Otherwise, A^^ —(do + di) is an odd 
number, so (59) is equivalent to 2d* > -f do + di — 1. 
In addition, Nb + do d- di owns the same parity as 
Nb — do — di, thus Wh + do + di — 1 is an even number. 
Therefore 2d* > Nb + do+ di. To sum up, (59) indicates 
2d* > A^h + do + di. Thus d* — do > A^b — d* + di. 
Moreover, to satisfy (30c), we should have A), — d* + 
di > rank{W}. So d* — do > rank{W}. However, 
rankjGiV} = rankjV} — do due to rankjV} — 
do < min{A^a, .^e}- Thus, rank{GiV} = d* — do > 
rank{W} > rank{H 2 W}, which contradicts (30b). 

In Case IV, d* = do + d 2 , where do = {Na — Ng)~^ and 
d 2 = min{s 4 , J}. Since the analysis is similar to Case 

III, so in the sequel we only give the skeleton on it. 

1) For the case of S4 < d* = do + S4, 

which, combined with the assumption d* > d*, gives 
d* > do + S4. Thus d* — do > S4. However, with the 
GSVD Transform of (H^, Gf^), s = dim{span(H 2 ) fl 
span(Gi)} = S 4 . To satisfy (30b), we should have 
rank{GiV} <5 = 54, which, together with the fact 
rank{V} — do < minlWo, Ng}, indicates d* — do < S4. 
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2) For the case of 54 > J, d* = do + 

which, together with the assumption d^ > d*, gives 
2d* > Nf, + do- Thus d* — dp > Nt, — dK However, 
to satisfy (30c), we should have iV;, — d* > rankjW}. 
Therefore, rankjV}—dp = d*—dp > rankjW}, which, 
together with the fact rankjGiV} = rankjV} — dp, 
contradicts (30b). 

Summarizing the above four cases, for any feasible points 
for the problem of (30), denoted by {V,W} and V G 
^NaxKa^ if Ka = d*, d* < d*. On the other hand, if 
Ka > d*, resorting to the singular value decomposition 
(SVD) of HiV, we can always find another feasible point 
{V', W'} for the problem of (30), such that V' G 
and = rank{V'} = rank{HiV'} = d*. As such, the 
assumption d* > d* also contradicts the feasibility conditions 
in (30). So d* < d*. 

In conclusion, for any feasible points for the problem of 
(30), we should have d* < d*. By definition, we arrive at 
d < d*. This completes the proof. 

Appendix G 
Proof of Corollary 2 

According to Table II, it is obvious that s.d.o.f. = 0 can 
only happen in the last case, where 

Na<Nb + Ne- Nj,Nb < N, <Ne+Nb-Na 
or Na<Ne,Nj <Nb. (60) 

Thus, to complete the proof, we only need to focus on the 
case of (60), where s.d.o.f. = min{s, ["^J} in which s = 
mm{Na,Ne} + mm{Nj,Ne} — mm{Na + Nj,Ne}. Note that 
the case Nb+Ne — Nj < Na < Ng, Nt < Nj < Ng+Nt — Na 
never happens, so (60) is equivalent to 

Na < Ng, Nj <Ng + Nb- Na, (61) 

which implies s = Ng + min{A(;, A’e} — min{A)j + Nj,Ng}. 
Therefore, 

1) for the case of Ng > Na + Nj, s = 0 thus s.d.o.f. = 0; 

2) for the case of Ng < Nj, s = Na thus s.d.o.f. = 
minjiVo, L^J}; 

3) for the case of Nj < Ng < Na + Nj, s = Na + Nj — Ng 
thus s.d.o.f. = min{Afo + Nj — Ng, 

In conclusion, when Nb = 1, s.d.o.f. = 0 if and only if 
Ng > Na + Nj — 1; when Nb > 1, s.d.o.f. = 0 if and only if 
Ng > Na + Nj. This completes the proof. 
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